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ABSTRACT 


Through the study of a carefully selected sample of isolated spiral galaxies, 
we have established that two important global physical quantities for tracing 
star forming activities, Lpj^i and have non-linear dependence on another 

commonly cited global quantity Lp- Furthermore we show that simple power 
law relations can effectively describe these non-linear relations for spiral galaxies 
spanning four orders of magnitude in FIR and Mh 2 and nearly three orders 
of magnitude in While the existence of non-linear dependence of Mh 2 
(assuming a constant CO-to-H 2 conversion) and Lpjp on optical luminosity 
{Lb) has been previously noted in the literature, an improper normalization 
of simple scaling hj Lp has been commonly used in many previous studies 
to claim enhanced molecular gas content and induced activities among tidally 
interacting and other types of galaxies. We remove these non-linear effects using 
the template relations derived from the isolated galaxy sample and conclude 
that strongly interacting galaxies do not have enhanced molecular gas content, 
contrary to previous claims. With these non-linear relations among Lp, Lpjp 
and Mh 2 properly taken into account, we conhrm again that the FIR emission 
and the star formation efficiency {Lpip/Mu^) are indeed enhanced by tidal 
interactions. Virgo galaxies show the same level of Mhj and LpiR as isolated 
galaxies. We do not find any evidence for enhanced star forming activity among 
barred galaxies. 


Subject headings: galaxies: fundamental parameters — galaxies: interactions 
galaxies: ISM — infrared: galaxies — radio lines: galaxies 
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1. Introduction 


The molecular gas (H 2 ) and far-infrared (FIR) contents are important quantities in 
analyzing the star formation activity in galaxies. They have been analyzed in the literature 
for different galaxy samples ([Young fc Scoville 1991| and references therein, ^ofue et al.| 
1993[ [Braine fc Combes 1993|, [Young et al. 1996| , Solomon et ah 19^) although largely 


among the IR luminous galaxies. These two quantities correlate linearly assuming a 
constant CO-to-H 2 conversion factor, and massive star formation within molecular clouds is 
sufficient to explain the observed trend in most cases. On the other hand, several authors 
have noted that these quantities do not correlate linearly with the total blue luminosity (Ls) 
in the sense that and Lfjr/Lb increase with increasing ([Young 1987| , [Young| 

et al. 19'89| , ^age fc Solomon 1989|) . Nevertheless, a simple normalization using Lb has 
been commonly used to claim enhanced activities among different groups of galaxies (e.g., 
Braine fc Combes 1993| , [Combes et al. 1994| , Pevereux fc Hameed 199'7[ ) while what they 
have actually shown is this residual dependence on the luminosity. In this paper, we first 
examine the nature of non-linearity involving Lb, and then re-examine the previous claims 
of enhanced molecular gas content and star forming activity based on the normalization 
by Lb- We show here that the enhanced molecular gas content among bright interacting 
galaxies reported in the literature is entirely due to the non-linearity effect, but the excess 
in the FIR emission is real and is largely due to real enhancement induced by environmental 
effects. 


We have compiled from the literature an extensive comparison sample of 207 galaxies in 
varying interaction classes and environment (isolated, weakly perturbed, strongly perturbed, 
Virgo cluster), representing a wide range of optical luminosity (10®'® Lq < Lb < 10^^'"^ Lq) 
and for which FIR, CO, and optical luminosity (as measured by the blue luminosities Lb) are 
available. The FIR luminosity is obtained as log{LFiB/ Lq) = log{FIR) + 2log{D) + 19.495 
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where D is distance in Mpc and FIR = 1.26 x 10 x (2.58/60 + -^loo) W m ^ ( |Helou et al. 


|1988|) . The 60 pm and 100 pm fluxes are obtained from the NED^ database. The molecular 
hydrogen mass is derived as Mh 2 = 4.82 aIcod% assuming a constant CO-to-H 2 
conversion factor of a = NH 2 /Ico = 3 x 10^° cm“^ (K km where Ico is the velocity 

integrated CO intensity in K km s“^ and is the half-power beam diameter in pc at the 
source distance ( [Sanders et al. 19911) . The use of a constant conversion factor is discussed 
further in § 2. Because many nearby galaxies typically subtend several arcmin in size and 
are much larger than the beam of the telescopes used, the comparison data consist mostly 
of the CO surveys using at least partial mapping. The optical luminosity is derived from 
the magnitude from the RC3 catalogue, corrected for galactic absorption (using the 
extinction value given by [Burstein fc Heiles 1984| , with the reddening law from [Savage ^ 
Mathis 19791) and internal extinction ([de Vaucouleurs et al. 19M[ RC3) using the redshift 
given in NED. We adopted a value for the Hubble constant of Hq = 75 kms“^ Mpc“^. 


2. Relationship among Lb, L^jb and Mh 2 for Isolated Galaxies 


In order to first understand the intrinsic non-linear nature of with respect to FIR 
luminosity and molecular gas content, we have assembled from the literature a reference 
sample of isolated galaxies representing a wide range of luminosity that extends from Lb 
= 10®'® to 10^°'® Lq. Our isolated galaxy sample consists of 68 objects from the distance 
limited survey of the Nearby Galaxies Catalog ([Tully 1988|) by [Sage (1993)[ and the class 
0 objects of the IRAS selected sample by [Solomon fc Sage (1988)[ . For six galaxies only 


^The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 
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upper limits are available on the CO luminosity. Morphological types range from Sa to Sd 
galaxies, each distributed along the entire luminosity range, and the types around Sc are 
the most abundant. This sample lacks completeness because the galaxies studied in the 
literature are frequently biased towards infrared luminous galaxies. On the other hand, our 
isolated galaxy sample is sufficient for this comparison study in the sense that the range of 
optical luminosity present in other comparison samples are basically included. 

In Fig. la, we plot Lb against Mhj for this sample. The data can be well described by 
a simple power law relation 

loglB = (0.57 ± 0.03)+ (4.9 ±0.6) (1) 

The coefficients are determined by a linear regression taking into account the upper limits 
using the Astronomy Survival Analysis package (ASURV^. We do not hnd any dependence 
on the morphological types. Sage k. Solomon (1989) report a slope of 0.53 ± 0.07 for a 
sample of isolated and weakly perturbed galaxies, and Young et al. (1989) hnd a slope 
of 0.72 ± 0.03 for a sample of 124 galaxies with CO data. Removing the 5 nearby dwarf 
galaxies with lowest luminosities (L^ < 7 x 10^ Lq) that clearly separate from the overall 
tendency, we hnd a slope of 0.54 for the Young et al. sample - which is consistent with our 
result. 

In Fig. lb we plot LpiR against Mhj, whose ratio is a measure of the star formation 
efficiency (SFE) (see [Young fc Scoville 199 1|) . This correlation has been widely considered 
in the literature and found to be nearly linear when a large range of luminosities are 

^Astronomy Survival Analysis (ASURV) Rev 1.2 package is a generalized statistical 
analysis package which implements the methods presented by Feigelson & Nelson (1985) 
and Isobe et al. (1986) and are described in detail in Isobe & Feigelson (1990) and La Valley 
et al. (1992) 
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considered, 

logLpiR = (0.90 ± 0.05) logMn^ + (1.5 ± 0.4) (2) 

with a similar dispersion as that of relation. 

These plots demonstrate that increases more slowly compared with the molecnlar 
gas content (for a constant CO-to-H 2 conversion factor) and far-infrared Inminosity, and 
more importantly, these non-linear dependencies are present even among these carefnlly 
selected isolated spirals. While the physical canse for the observed non-linear relations 
is nncertain, it is snggestive that the physical mechanism responsible for the FIR and 
CO emission may be closely related, bnt they are distinct in natnre from optical 

emission. For the ease of the snbseqnent discnssions, we have parameterized the above 
empirical non-linear relations in terms of Mh 2 , bnt we do not claim that Mh 2 (or any of the 
three physical qnantities compared here) is the fnndamental driver. 

If one postniates that the root of the non-linear relations shown in Eqs. 1 & 2 lies in 
the FIR and CO emission, then an obvions possible canse is the metallicity dependence on 
the Inminosity. In fact, a positive correlation between metallicity and absolnte magnitnde 
has been obtained by Roberts & Haynes (1994) for a large sample of spiral galaxies, 
and Arimoto et al. (1996) have proposed an empirical dependence between CO-to-H 2 
conversion and absolnte magnitnde based on their observational data and the CO emission 
model of Sakamoto (1996). Indeed there are clear indications that the CO intensity may 
not be a good tracer of the molecnlar content among metal-poor galaxies (|Cohen et al.| 
19881, [Rnbio et al. 1991]) . Also, the non-linear dependence in Eq. 1 is nearly reprodnced by 
the metallicity dependence derived by Roberts & Haynes if the CO emission is assnmed to 
depend linearly on the metallicity. However, there are observational and theoretical reasons 
to believe that snch a metallicity dependence cannot work. For example, the observational 
evidence for metallicity dependence offered by Arimoto et al. becomes marginal if the low 
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spatial resolution data on M51 and LMC are removed. Also, full chemical and radiative 
transfer modeling of CO emission by Wolfire, Hollenbach, & Tielens (1993) and Maloney 
& Wolfire (1997) snggests that CO emission is a robust tracer of Mh 2 against metallicity 
variation - no metallicity dependence is expected nntil the metal abnndance becomes low 
enongh to make the CO transition optically thin. 

Alternatively if the non-linear relations derived above are rooted in L^, then optical 
extinction and scattering model may offer an explanation. A naive expectation is that 
larger extinction in more massive galaxies with larger dust column density would result in 
a slower increase in Lb with size (e.g., Young et ah 1989). Bnat & Xn (1996) also hnd a 
good correlation between the optical depth in B and the ratio of FIR to blue luminosity. 
Therefore, it is natnral to suspect that the root of the observed non-linearity may lie with 
Lb rather than with Lf/j? or 

In a simple geometrical extinction model in which optically thick tracers snch as 
blue light represent surface area and optically thin tracers snch as FIR and integrated 
CO intensity represent volnme elements, a non-linear dependence is natnrally expected 
between these two classes of tracers - i.e.. Lb oc for a spherical nniform clond. 

In reality, the structure of the ISM is more complex and clumpy or hlamentary, and the 
effective optical depth is signihcantly reduced compared with a uniform screen model by 
extinction geometry and efficient scattering (Witt et al. 1992, 1996). Nevertheless, it is 
remarkable that the non-linear relation between Lb and Mh 2 is observed with nearly the 
same power law dependence as predicted by this simple model, and one may infer that the 
basic assumptions of the model stated above must hold at the clump-interclump size scales. 
We note that the predicted exponent of 2/3 is independent of the clond shape as long as 
the clond structure remains self-similar at different size scales. 

The extinction model also correctly predicts other optical depth dependent trends. For 



example, the exponent of the correlation between Lb and Lpm is expected to be steeper 
than that of Lb-Mh 2 because additional surface element provided by HI among the low 
luminosity galaxies is partly compensated in FIR by the dust emission associated with the 
HI gas. A formal £t for the Lb-Lfir relation in our isolated galaxy sample gives 

IoqLb = (0.65 ± 0.09) loglpm + (3.9 ± 0.9) (3) 

and the slope of the correlation is indeed larger than that of Lb-Mh^ (0.57 ± 0.03). We 
predict that the optical luminosity at longer wavelengths {Lr or Lk) would produce a more 
linear relationship with Mr^ than Lb- 


3. Re-evaluation of the Tidal and Environmental Effects 

Regardless of their physical causes, two global physical quantities LpiR and Mh 2 show 
clearly non-linear dependence with respect to among a large sample of isolated galaxies. 
The relations are well described by a simple power law spanning four orders of magnitude 
in Lfir and Mhj and nearly three orders of magnitude in (see §2). Therefore, any 
comparison of Lfir and Mh 2 normalized simply hy Lr produce a misleading result. This 
is particularly true for the studies of bright interacting galaxies found in the literature 
where the samples are often selected based on their high luminosity. Here we re-evaluate 
the impact of tidal and environmental effects on individual galaxies by examining their 
global properties using the empirical relations derived from our isolated galaxy sample as 
templates in order to discriminate the size-dependent effects from the environmentally 
induced effects. Deviations from the template power laws (Eqs. 1 & 2) are evaluated 
below for galaxies suffering different degrees of interactions. We also examine whether the 
presence of a bar correlates with any enhancement in the star formation efficiency. While 
each subsamples analyzed here represent different range of luminosity with different mean 
luminosity, almost the entire range of optical luminosity considered are represented in our 
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isolated galaxy sample, and the use of these template relations are justihed. 

Both the FIR and CO data available in the literature are severely sensitivity limited, 
and a special care has to be made to incorporate the upper limits properly in the analysis. 
Standard statistical tests such as Kolmogorov-Smirnov (KS) or Mann-Whitney cannot 
be rigorously applied to these heavily censored data, and we instead use equivalent tests 
provided by ASURV. In particular to quantify the statistical probability of the compared 
samples representing the same distribution function, we use the logrank and Peto-Prentice 
generalized Wilcoxon test which are known to be the most robust against the censoring 
pattern. 


3.1. Interacting and Cluster Galaxies 


The entire sample studied here consists of 139 galaxies, divided into three subclasses: 
weakly perturbed (WP), strongly perturbed (SP), and Virgo Cluster (VC) galaxies. The 
WP sample has 43 galaxies including interaction class 1, 2 & 3 of |Solomon fc Sage (1988) 


and class 2 objects from the luminous IRAS sample by Sanders et al. (1991). The SP 
sample has 38 galaxies including interaction class 4 of [Solomon fc Sage (1988)| , interaction 
class 3 & 4 of Sanders et al. (1991), and closely interacting pairs from Combes et al. (1994). 
The dehnitions of “weakly perturbed” and “strongly perturbed” are given in the references 
listed above, SP galaxies being those identihed as the hnal stages of mergers. In addition, 
we have constructed a sample made of 58 member galaxies in Virgo cluster, including both 
bright (Kenney & Young 1988ab) and faint spirals (Boselli et al. 1995). In the SP and VC 
sample, 3 and 18 galaxies have only upper limits in CO, respectively. 


In Fig. 2, we plot Mh 2 against Lb for the VC, WP, and SP samples along with a solid 
line representing the power law derived from our isolated galaxies sample (Eq. 1). We have 
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analyzed the residuals of each subsample with respect to the isolated galaxies template, 
and none of the three subsamples can be distinguished from the isolated galaxy template. 
In Table 1, we summarize the median values of the residuals and the semi-interquartile 
distances which measure the dispersion of the distribution. The deviation of the residuals 
from zero value are statistically negligible, and this fact indicates that neither the optical 
luminosity nor the molecular gas content as measured in CO are signihcantly affected by 
tidal interactions. This contradicts the previous reports of enhanced molecular gas content 
among strongly interacting pairs (Braine & Combes 1993, Combes et al 1994). Virgo spirals 
show a gas content similar to isolated galaxies as has been previously suggested (see Boselli 
1994 and references therein). 

The relation between LpiR and Mhj is plotted in Fig. 3 for the three subsamples 
along with the power law relation obtained from the isolated galaxy sample. While the 
Virgo spirals and WP sample galaxies follow the same Lfir-Mhj relation as the isolated 
galaxy sample, the SP sample galaxies show a clear deviation - enhanced FIR emission for 
given molecular gas content, already known in the literature as increase “star formation 
efficiency” (Lpiji/Mh 2 ratio - e.g. Sanders et al. 1991). The FIR luminosity of the SP 
sample galaxies lies mostly within the range covered by the isolated galaxy sample, but it 
extends to higher FIR luminosity and has a higher mean FIR luminosity than the isolated 
sample, requiring some extrapolation of Eqs. 1 & 2. We note however that the ranges 
of Mh 2 and Lr covered by the SP and isolated comparison sample are the same, clearly 
indicating that the FIR emission is enhanced among the SP sample galaxies. We have 
improved the understanding of Lpjp/Mh^ ratio by proper accounting the strong non-linear 
luminosity dependence while previous studies of Lpjp/Mp^ ratio were performed with a 
simple assumption of no luminosity dependence and despite claims of enhanced molecular 
gas content in some cases. 
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Since we have shown that H 2 content is independent of tidal disruption or environment, 
any deviations from the template will indicate changes in the FIR emission. In Table 2 
we summarize the statistics and probability of association in distribution for 

the VC, WP, and SP samples with respect to the isolated galaxy sample. A zero value 
for the probability indicates that the compared distributions are statistically different 
from the isolated galaxy sample. Also included in Table 2 are the median and the 
semi-interquartile distance of the log[LFiR/MH 2 ] for all the samples as obtained from the 
cumulative distribution functions. All performed statistical tests indicate large differences 
in the distribution of the SP galaxies with respect to the isolated galaxies, as it was already 
clearly shown in Fig. 3. The VC galaxies show a Lpjji-Mh 2 distribution very similar to the 
isolated galaxies as indicated by the similarity of their median values and by high statistical 
likelihood. 


An evident result is that the star formation efficiency in the SP sample is enhanced 
on average by a factor 3.5 with respect to the isolated galaxy sample, consistent with 
the results in the literature (ganders et ah 1986i Solomon & Sage 1988, [Young fc Scoville| 
1991 ), even when various non-linear dependences are properly accounted (the probability 


of the SP sample having the same distribution as the isolated sample is zero). Our new 
analysis also removes any uncertainty regarding variations in Mh^ as function of interaction 
strength (thus affecting Lpjr/Mh 2 ratio) since no such variation is found. The weakly 
interacting systems show some marginal excess in SFE (30%, well within one sigma). This 
can be understood since observational separation between weak and strong interactions is 
not always straightforward, and some contamination exists between the two samples. 
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3.2. Barred and unbarred galaxies 


Out of 196 galaxies in our sample with morphological classihcation listed in NED, 117 
are classihed as barred (38 isolated, 28 WP, 13 SP, and 38 VC) and 79 as unbarred (30 
isolated, 15 WP, 14 SP, and 20 VC). The same analysis is now performed for the barred 
versus unbarred galaxies in order to test the previous claims that enhance star forming 
activity is associated with barred spirals (e.g., [Young 1992|) . 


In Fig. 4a, we plot Lb against Mh 2 for the barred and unbarred galaxies, and no clear 
difference is seen between them, indicating that the total molecular mass is independent of 
their barred morphology. In fact Lb distributions of barred and unbarred galaxies are found 
indistinguishable by different statistical tests applied, and the same conclusion is drawn 
for their Mhj distributions also (< Ioq^Lb) >bar = 10.231, < Ioq^Lb) >unbar = 10.110; 

< log[MH 2 ] >bar = 9.195, < log[MH 2 ] >unbar = 9.110). 


We have also explored for any differences in distribution as shown in Fig. 

4b. Again the data overlap, and the FIR emission of the barred and unbarred galaxies 
is found indistinguishable by the statistical tests performed (< log^LpjB) >i,ar = 9.691, 

< log{LpiB) >unbar = 9.687). Thus the star formation efficiency seems independent of the 
barred morphology, as shown in the histogram of the ratio Lpm/M h 2 in Fig- 5. The median 
values that are respectively [LFjR/MH 2 ]unbarred = 0.69 {a = 1.59) and [LFiR/MH 2 ]barred = 
0.76 (cr = 1.58). Since bar morphology and starburst activity are sometimes associated with 
tidal interactions, this result may indicate that the lifetime of a stellar bar, if transient, is 
much longer than the time scale for starburst. 



4. Conclusions 


Through the study of a carefully selected sample of isolated spiral galaxies, we have 
clearly established that the two important global physical quantities for tracing star forming 
activities, Lpjji and have non-linear dependence on another commonly cited global 

quantity Lp- Furthermore we show that simple power law relations can effectively describe 
these non-linear relations for spiral galaxies spanning four orders of magnitude in FIR and 
Mh 2 and nearly three orders of magnitude in These non-linear relations are critically 
important in comparing different samples of galaxies since a simple normalization and the 
resulting residual dependence imply that samples consisting of the more luminous galaxies 
would show intrinsically larger ratios. Using the non-linear template relations obtained 
from the isolated galaxy sample, we have re-analyzed previous claims of enhanced molecular 
gas content and star formation through tidal interactions. We have also examined the effect 
of cluster environment and presence of stellar bars on the H 2 content and any enhancement 
of star formation. 

In the frame dehned by the optical luminosity, molecular gas content, and far infrared 
luminosity, we found that gravitational interaction affects only the FIR emission. We found 
no evidence for enhanced molecular gas content among the interacting galaxies, and we 
conclude that the previous such claims are due to the use of size normalization using Lp 
without properly accounting of the non-linear dependence. We do conhrm that enhanced 
FIR emission (star formation) by tidal interaction is real. We found no evidence for 
enhanced star formation among barred galaxies in comparison with the unbarred galaxies. 

The authors thank N. Scoville, P. Solomon and J. Sulentic for useful discussions. JP, 
AO and LV-M are partially supported by DGICYT (Spain) Grant PB93-0159 and Junta de 
Andalucia (Spain). MSY is supported by the NRAO Jansky Research Fellowship. 
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Fig. 1.— The dependence of Lb and Lpm on Mhj for the isolated galaxy sample. The solid 
lines correspond to the htted power law given in Eq. 1. 

Fig. 2.— The dependence of Lb on Mhj for the three samples of interacting galaxies. The 
plotted line corresponds to the template power law derived from the isolated galaxies. 

Fig. 3.— The dependence of LpiR on Mhj for the interacting galaxy sample. 

Fig. 4.— The dependence of Lb and Lpm on Mhj for all galaxies in our sample, separated 
by barred and unbarred morphology. 

Fig. 5.— Histogram of the star formation efficiency measured by the ratio [Lfir/Mh. 2 \ for 
our full sample, divided into barred and unbarred galaxies. 
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Sample 

Median 

Q“ 

Isolated 

0.043 

0.17 

W. Per. 

0.204 

0.21 

S. Per. 

0.101 

0.27 

Virgo 

0.197 

0.14 


Table 1: Statistical parameters of the residuals relative to the Mh 2 ^Lb template obtained 
from isolated galaxies. 


“ Semi-interquartile distance. For a normal distribution a = 3/2 Q. 




Isolated 

Virgo 

WP 

SP 

logrank statistic 

— 

1.103 

2.193 

5.769 

Probability 

— 

0.270 

0.028 

0.000 

Peto-Prentice 

— 

0.191 

1.889 

6.457 

Probability 

— 

0.848 

0.059 

0.000 

Median 

0.633 

0.583 

0.737 

1.171 

Semi-interquartile 

0.287 

0.404 

0.279 

0.184 


Table 2: Statistical parameters of log[LFiR/MH 2 ] distribution for the studied samples as 
compared with the isolated galaxies template (from Astronomy Survival Analysis package). 























